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Abstract 

(~-^ ■ Wigner phase space quasi-probabihty distribution function is a Fourier transform related to 

. a given quantum mechanical wave function. It is shown that for the wave functions of type 

, '^{q) = e~"'^ </'('?), the Wigner function can be defined in terms of differential operators acting 

^ ' on a given function, independently from the integral formula which appears in the standard 

QJ • definition. Gaussian wave packet, harmonic and singular oscillators are given as the examples. 

\0 ; PACS: 03.65.-W, 03.65.Sq, 05.30.-d 
> 

oo 

o 

^ : 1 Introduction 
o 



X 



In a study of quantum corrections to classical statistical mechanics, Wigner constructed the function 

W{q,p) = —j riq-^y)HQ + ^y)e-'''dy. (1) 

^ I that now called by his name [1]. The goal was to replace the quantum mechanical wave function 
ip{q) with a probability distribution in phase space that providing a particlelike description of the 
underlying wave propagation. Hence the Wigner function (WF) exhibits a representation of quantum 
states in phase space. Although in quantum mechanics due to the uncertainty principle, a certain 
point in phase space does not make sense, projection of the WF onto the q and p planes generates the 
relevant marginal distributions associated with that quantum state, where possible negative values 
correspond to non-classical situations |2]. 

WF has been much studied theoretically and experimentally since its introduction, not only in 
the context of quantum mechanics, but also in various branches of physics and related sciences. 

The WF given by ([1]) can be obtained easily if the integral allow to get an answer analytically. 
In most cases, despite the existence of wave function tp{q), the integral becomes quite cumbersome 
and sometimes it is impossible to handle it. 

The aim of the present paper is to present an alternative definition of the WF for the wave 
functions of a certain type, that is a derivational approach rather than the original integral based one. 
This new definition gives, at least partly, calculational tricks to overcome the difficulties mentioned 
above. 
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The approach considers a class of wave functions given in the form 

m = e-"^V(g) (a ^ 0), (2) 

where a is a positive constant, and is based on the series expansion of (f){q) in term of Hermite 
polynomials. In order to establish the convergence of the series we need the following theorem 
If 0(g) is a piecewise smooth real function in every finite subinterval [—6, b] and if the integral 

/oo 
e-'^'<P\q)dq (3) 
-oo 

is finite, then the series 

oo 

= '^nHn{q), -oo < g < oo, (4) 

n=0 

where the coefficients C„ can be computed by the orthogonality properties of the Hermite polyno- 
mials, converges to 0(g) at every continuity point of 0(g) . It is also assumed that in the case of 
complex 0(g) the consequent of the theorem is still valid i.e., (j){q) admits a uniformly convergent 
expansion on the basis of Hermite polynomials. 

With these assumptions it is possible to replace the integral by the series and finally to define 
the WF as the following 

W{q,p) = e-^'^^V* (g - ^S,) (g + e'^V^a.^, (5) 

which is the main result of this paper. Note that the assumptions ensure that 0(g) is smooth enough 
to be evaluated over the differential operators. It is easy to show that the operators and 0* in are 
commutative and their product is real. A useful aspect of this definition is that if 0(g) contains a finite 
polynomial alone that may be an orthogonal polynomial such as Hermite or (associated) Laguerre 
polynomial, then obtaining the WFs for the individual states requires only simple differentiations 
without needing integral tables or computers. 



2 Derivation 

With the choice of ([2D the WF ([I]) amounts to 

W{q,p) = —e-^'^'^' / dy <P*{q- -y) <P{q+ -y) e-'-^^y'-^y^ dy. (6) 



For our purpose, the Hermite polynomials takes the form 
By the substitution of (jl]) and ([7]) into (JH]) we get 

w{q,p) = -e-'^'^' y c:c^ y ^- o-^-y ^> ^- o-^- 

^^'^^ 27r « r\ (n - 2r)\ ^ s\ (m - 2s)\ 

n,m=0 r=0 ^ ' s=0 ^ ' 

{q - ^yr-'''{q + ^yT''' e-'-^^y'-^y^ dy. (8) 
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If one uses the Binomial expansion for the terms in the integral in ([8]), the WF takes the form 

n,m=0 r=0 ^ ' s=0 ^ ^ 

/•oo 2 

X / 7/"+'^ e-^'^^/^-'^P dy, (9) 



x(p) 

where the integral can be adopted to the integral p] 

[f e-^+--cit = e-^^if„(x). (10) 

Thus is straightforward; 

X{P) = (-O'v^^ (^^j e-^^/2'^^^i/,(p/^V^), (11) 

where A; = a + /?. If we use the Rodrigues representation of the Hermite polynomials 

H,iu) = i-lfe^'d':e-^\ (12) 

we get x{p) ^ 

_ 

Thus the WF yields 



x{v) = ^^W^'e-^''-^^\ (13) 



X 

/3=0 



.:^o " "^^H(n-2r)! ^.!(m-2.)! 
""^^M g'"-'^-^(|)^(29p)^e-^''/2«?^'. (14) 



The last two sums in ([H]) stand for (g - ySp)" and (g + ySp)"" respectively. With this 
arrangement it is obtained that 



^ oo [n/2] / I -fc 

* n=0 r=0 



oo [m/2] / i\s I •»- 

X Ec-Epfe^pfe+ff.)!-"^-'"'"""- (15) 

m=0 s=0 ■ ^ 

Thus, with the help of (jl]) and ([7]), the Wigner function can finally be compacted as in the equation 
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3 Applications 

3.1 Gaussian wave packet 

The wave function for a system represented initially by a Gaussian wave packet is given as 



ijj{q) = C exp 



4(Ag)2 _ 



(16) 



where C = l/[27r(Ag)2]^/^, Ag is the width of the packet centered at go and po is its average 
momentum. According to ([2]) 



0(g) = C exp 



ql 



4(Ag)^ 



(17) 



and a = l/[4(Ag)2]. Thus 



^*\q- -dp 



~2 



q + y (9p ) = exp 



ql 



;i8) 



which is obviously a translation operator which converts p to p — Pq. Therefore (j5]) and the fact that 
AgAp = h/2 yield the WF as 



1 


{q-qof 




{P-Pof 


— exp 
nn 


2(Ag)2 _ 


exp 


2(Ap)2 _ 



that confirms the correct result. 



(19) 



3.2 Harmonic oscillator 

The wave function corresponding to the Hamiltonian H = p'^/2 + g^/2 is given by 

Mq) = Cne~'^"/^''H^{q/Vh), 



where C„ = [1/ (nh)Y^'^/ V2"-n\ . By virtue of the equation [1] 



H„ 



1 f '^f, 



1 

71 



q+i^d. 



= (-1)" 2" n! L„[(2g2 + 2p^)/h] e"^'/'' 
which is obtained by iteration, gives the well known WF 

-1^ 



Wniq,p) 



TT 



h 



e-'^'/^Lnim/h), 



(20) 



(21) 



(22) 



where H = p'^/2 + g^/2 and L„ denotes the Laguerre polynomial of order n. 
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3.3 Singular oscillator 



The system is described by the Hamiltonian H = /2 + (f /2 + g'^ / (f with g is a, real constant. The 
normahzed eigenf unctions of the Hamiltonian are 

i,n{q)=Cnq''K-"\q'/h)e-'^'"\ (23) 

where = n\/\r{n+a+l/2) K^^^], T is the Gamma function, L"^^^^ is the associated Laguerre poly- 
nomial and a = 1/2 + (1/4 + 2(7^)^/^ l^j. The operator equation version of the WF is straightforward 
with the help of the argument presented up to now; 

... / r^2^ 

Wn{q,p) 



idp 



idr. 



(24) 



For an arbitrary a and n, integral tables and the approach discussed here fail to get an exact WF 
of this system. Especially for the fractional values of a, even determining the ground state carries 
big difficulty since the action of the operator (g^ + 9^/4)" on exp{—p'^/h) is unknown. (At least, I 
am not aware of it). But for some special values of a and n, WFs for this system can be obtained 
explicitly [4J. 

As a limit of the applicability of the method, consider the one dimensional anyon system [6], 
where the wave function 



(25) 



corresponds to the Hamiltonian H = p^ /2 — 1/q — ^/(f- Obviously the matter arises from the term 
g^/^ since the smoothness condition of 0(g) is violated. 
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